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The largely unexplored decameter ra- 
dio band (10-30 MHz) provides a unique 
window for studying a range of astro- 
nomical topics, such as auroral emission 
from exoplanets, inefficient cosmic ray 
acceleration mechanisms and fossil radio 
plasma. The scarcity of low-frequency 
studies is mainly due to severe iono- 
spheric corruptions. Here we present 
a calibration strategy to correct for the 
ionosphere in the decameter band. We 
apply this to an observation from the 
Low Frequency Array (LOFAR) between 
16-30 MHz. The resulting image cov- 
ers 330 square degrees of sky at a res- 
olution of 45”, reaching a sensitivity of 
12mJy beam~!, improving by an order 
of magnitude in terms of sensitivity and 
resolution compared to previous decame- 
ter observations. Residual ionospheric 
effects cause additional blurring between 
60-100”. In the surveyed region, we 
have identified four fossil plasma sources. 
These sources likely harbor rejuvenated 
radio plasma from past active galactic 
nuclei outbursts. Three are situated near 
the center of low-mass galaxy clusters. 
Notably, two of these sources display the 
steepest radio spectral index among all 
the sources detected at 23 MHz. This 
indicates that fossil plasma sources con- 
stitute the primary population of steep- 


spectrum sources at these frequencies, 
emphasising the large discovery potential 
of ground-based decameter observations. 

The first detection of extraterrestrial radio 
waves was made in the decameter band by Karl 
G. Jansky [i]. After these initial observations, 
most studies of the radio sky were carried out at 
higher frequencies due to greater observational 
and technical advantages, such as increased res- 
olution, lower sky temperature, and improved 
ability to be calibrated. Therefore, the decame- 
ter band remains a relatively unexplored spec- 
tral window. There has recently been a resur- 
gence of interest in the decameter band because 
it is a unique window for studying several astro- 
physical sources and radiation mechanisms that 
are inaccessible at higher frequencies. Some 
classes of radio sources are known to emit pri- 
marily at low frequencies with a sharp spec- 
tral cutoff. Others have ultra-steep radio spec- 
tra throughout the radio frequencies at which 
they have been observed. Furthermore, sev- 
eral absorption processes which suppress low- 
frequency emission become more dominant in 
the decameter band [e.g., free-free absorption, 
synchrotron self-absorption; [3]. Decame- 
ter observations can provide critical constraints 
on the underlying emission and absorption pro- 
cesses. 

Coherent radio emission from stellar sys- 
tems can be bright in the decameter band [4]. 
This coherent emission is produced by plasma 


or electron cyclotron maser instability (ECMI) 
and is often generated in the coronae of stars. 
ECMI emission can also result from the in- 
teraction of a stellar or planetary companion 
with the magnetosphere of the primary body 
[4]. These auroral emissions are ubiquitous 
for planets in our Solar System which harbor 
strong magnetic fields, including the Earth. In 
the case of Jupiter, one prominent pathway 
by which ECMI emission can be induced is 
the interaction of Io with Jupiter’s magneto- 
sphere [5]. The radio spectrum from ECMI 
emission has a strong cutoff in the spectrum. 
For example, Jupiter is four orders of mag- 
nitude brighter in the decameter wavelength 
band than at frequencies above 30 MHz [6], 
which makes decameter observations uniquely 
capable for observing similar systems (e.g. in 
star-planet interaction systems [7]). Pushing 
the low-frequency limit is important to enable 
the detection of auroral emission from exoplan- 
ets with weaker magnetic fields since the cutoff 
frequency is directly proportional to the mag- 
netic field strength of the emitting body. The 
decameter radio band is also sensitive to the 
lower energy synchrotron-emitting cosmic ray 
electrons that trace long timescales and ineffi- 
cient acceleration mechanisms (such as turbu- 
lent re-acceleration [8]). For example, the lobes 
of radio galaxies will steepen with age when the 
active galactic nucleus (AGN) ceases its activ- 
ity, due to synchrotron and Inverse Compton 
losses. When this old (~ 10’~® yr) AGN radio 
plasma is located in the hot plasma that perme- 
ates galaxy clusters, it can be re-energised by 
adiabatic compression or other re-acceleration 
mechanisms [9] [0] [Ii]. This produces sources 
with complex morphologies that almost exclu- 
sively emit at low frequencies [I2]. Decameter 
observations thus provide a unique view of the 
complex life cycle of cosmic rays in galaxy clus- 
ters. 

While there are several important astronom- 
ical phenomena that can only be probed in 
the decameter band, there are three problems 
that must be overcome when observing at such 
frequencies. First, because of the long wave- 
lengths involved, the resolution that can be ob- 
tained is inherently limited. Achieving subar- 
cminute resolution requires baselines of the or- 
der of 100km. Second, the decameter band of- 
ten suffers from strong radio frequency inter- 
ference (RFI, e.g. [14]), particularly from 


shortwave radio communications. Third, the 
ionosphere significantly influences the propaga- 
tion of radio waves, which severely complicates 
the necessary high-resolution imaging. Vari- 
ations in the ionosphere’s total free electron 
content (TEC) cause variable phase delays in 
the arriving wavefront, in turn blurring radio 
images. In addition, the interactions of these 
electrons with magnetic fields in the ionosphere 
cause differential Faraday rotation. Both ef- 
fects are also time (on timescales of ~ seconds) 
and direction (on scales of a degree on the sky) 
dependent. The phase delays from TEC differ- 
ences scale linearly with the wavelength of in- 
cident radio waves, and Faraday rotation scales 
with the wavelength squared [e.g., [5]. There- 
fore, when observing at longer wavelengths, 
these effects become considerably more severe, 
ultimately becoming uncorrectable when ap- 
proaching the plasma frequency of the iono- 
sphere at around 3-10 MHz [e.g.[I6]. Given the 
above challenges, space-based radio arrays at 
frequencies below 30 MHz have been considered 
[e.g., [17]. Only a few decameter ground-based 
studies have been carried out that were sensi- 
tive enough to detect more than several dozen 
sources [e.g., [18] [9]. Surveys covering a sig- 
nificant part of the sky have been performed 
by the Dominion Radio Astrophysical Observa- 
tory (22 MHz, ref. [20]) and the Ukrainian T- 
Shaped Radio Telescope Second Modification 
(UTR-2, 10-25 MHz ref.[19]). The resolutions 
of these surveys range from approximately half 
to a few degrees, with noise levels of a few Jan- 
sky per beam or higher. Slightly outside the 
decameter band, the Eighth Cambridge Cata- 
logue of Radio Sources (8C, ref. [27]) at 38 MHz 
cataloged the sky above a declination of +60° 
at a resolution of approximately 4.5’. This sur- 
vey achieved a noise level of 0.3 Jy beam™ t. 
Recently, significant progress in the calibra- 
tion of low-frequency radio observations has 
been made with the Low-Frequency Array (LO- 
FAR, [22]), a pan-European radio array sensi- 
tive to frequencies between 10 and 240 MHz. 
The core of LOFAR comprises baselines up 
to 120km in the Netherlands. Previous stud- 
ies with LOFAR consist of, among others, the 
LOFAR, Two-meter Sky Survey (LoTSS, 120- 
168 MHz, ref. [23] (24]) and the LOFAR LBA 
Sky Survey (LoLSS, 42-68 MHz, ref. [25] [26]). 
The sensitivity, resolution, and frequency cov- 
erage of these observations, along with previous 
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Figure 1: Comparison of the RMS noise and central frequency of previous low-frequency obser- 
vations. Highlighted in blue is the wavelength range covered by the data described in this work. 
A Jupiter spectrum (extracted from Zarka et al. (2007) [54|) scaled to match the sensitivity of 
this work, is included to illustrate the importance of the decameter band for detecting electron- 
cyclotron maser emission. The size of the circles represents the resolution of the survey. For the 
LOFAR surveys (LoTSS, LoLSS), the frequency ranges are indicated with grey-shaded bands. 
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low-frequency work, are shown in 

In this work, we demonstrate the feasibility 
of sensitive subarcminute resolution decame- 
ter band observations using a calibration strat- 
egy to correct for the severe ionospheric image 
blurring that occurs below 30 MHz. We ap- 
ply this to a 5hr dual-beam LOFAR observa- 
tion between 16 and 30 MHz, going three times 
lower in frequency than LoLSS. This calibration 
strategy allows us to create a sub-arcminute 
resolution image of the decameter sky for the 
first time (see [Figure 2). To verify that the 
ionospheric conditions during our observations 
were typical, we analyzed observations of the 
primary calibrator over several nights. These 
demonstrated that the ionospheric conditions 
present during our observations were typical, 
based on the S4 index (see Methods). 

Our observing strategy consisted of simulta- 
neously observing a bright primary calibrator 
(3C 196) and the target fields. By scheduling 
the observation after midnight, we minimized 
RFI caused by the internal reflection of terres- 
trial RFI by the ionosphere, which is signifi- 
cantly worse during the day, as ionizing radia- 
tion from the Sun increases the column density 
of ions in the ionosphere. The remaining RFT is 
excised at high time (1 second) and frequency 
(3.05 kHz) resolution. The overall amount of 
data lost due to RFI was less than 10% (see 
Figure 3). 

By adapting the method of van Weeren et al. 
(2016) [27] and de Gasperin et al. (2021) [25], 
using individual facets calibrated for direction- 
dependent effects, we show that this approach 
can be even successful in the decameter radio 
band. First, instrumental effects are corrected 
by using 3C 196 as a calibrator source. Next, 
differential Faraday rotation and phase correc- 
tions for first and third-order ionospheric ef- 
fects were obtained for the whole field of view, 
after which the full field was imaged. This pre- 
liminary image reveals about 30 bright sources, 
although the majority of sources are severely 
blurred and smeared out, due to direction- 
dependent effects of the ionosphere. Follow- 
ing this step, we split up our field of view into 
several smaller “facets” and self-calibrated each 
facet individually, against a starting model ob- 
tained from the 151 MHz TGSS ADR1 [28] 
survey. This yields an improved image and 
model of the sky, partly corrected for direction- 
dependent effects. We repeated this step three 


more times, increasing the number of facets and 
iteratively constructing a better model of the 
sky, which in turn allowed us to calibrate on 
facets containing less flux. Based on the activ- 
ity of the ionosphere, 25 facets were calibrated 
and used in order to create a sharp image of 
the full field of view. 

The resulting image has a resolution of 45” 
at a central frequency of 23 MHz, and reaches 
a central noise level of 12mJy beam~', which 
is more than three orders of magnitude in im- 
provement compared to previous work below 
30 MHz {e.g. [19] [20]. The full region mapped 
is presented in The image covers 
an area of 330 square degrees, consisting of 
the sky area within 7 degrees of the pointing 
center of each observation. This corresponds 
to 0.6 times the full width at half maximum 
(FWHM) of the primary beam of the used ar- 
ray at 23 MHz. Beyond this radius, the image 
quality rapidly decreases due to the sparsity of 
the direction-dependent corrections. The im- 
provement due to our facet-based, direction- 
dependent approach can be seen in the cutout 
shown in [Figure 2| While the image quality is 
greatly improved, some artifacts remain around 
the brighter sources due to imperfections in the 
ionospheric corrections. These imperfections 
can be modeled with a Gaussian blurring of 
60” for brighter sources, up to 100” for fainter 
sources. 

A source catalog was constructed from the 
image using PyBDSF [29], with flux densities 
measured in apertures having radii of 2 arcmin. 
This ensures that the measured flux densities 
are accurate regardless of the shape of the point 
spread function. The source catalog can be 
found online, or with the paper as Supple- 
mentary Data 1. We detected a total of 2863 
sources within the 330 square degrees coverage 
of the image and a signal-to-noise ratio of at 
least 7, defined by the ratio between the flux 
density of a source divided by the local RMS. 
This corresponds to a source density of 8.7 
sources per square degree and an improvement 
of two orders of magnitude compared with the 
UTR-2 and DRAO 22 MHz surveys, the bench- 
mark surveys at these frequencies. Our mea- 
sured source positions are consistent to within 
< 10arcsec with those from TGSS and there 
appears to be no systematic flux-density scale 
offset compared with previous source catalogs 
(see Methods). 
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Figure 2: Overview of our observed region. Left: Full image of our observed region at a central 
frequency of 23 MHz. The image is clipped at a field of view of 7 degrees, which corresponds with 
the FWHM of the LOFAR beam at these frequencies. The black dashed rectangle corresponds 
to the region highlighted in the three panels on the right. Right: Comparison between direction- 
independent calibration (top) and the facet-based direction-dependent calibration (middle) of 
our data. The prominent ’spikes’ in the direction-independent image are caused by imperfect 
corrections of the ionosphere, as the ionosphere changes rapidly across the field of view. With 
direction-dependent calibration, the full field is divided into smaller facets that are individually 
corrected for the ionosphere, thus mitigating most of these spikes. As a comparison, the 8C 
image of the same area is given (bottom). 
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Figure 3: Fraction of LOFAR data flagged due 
to RFI as a function of frequency. The top 
image shows the observation used in this pa- 
per and the bottom image shows an observation 
taken in the early evening local time. There was 
significantly more RFI during the latter obser- 
vation, by about a factor of 7. Similar bad RFI 
conditions are experienced in the morning and 
during the day. 


The source catalog provided in this paper 
is used for constructing Euclidean normalized 
source counts in Our source 
counts show similar behavior to source counts 
obtained at higher frequencies, with a grad- 
ual decrease in the counts towards lower in- 
tegrated flux densities. From the figure, it is 
clear that the catalog is reasonably complete 
down to a flux density of around 200 mJy, be- 
low which the correction factors for incomplete- 
ness become significant. With an assumed blur- 
ring of around a factor 3 between the peak flux 
density and integrated flux density (see meth- 
ods), this corresponds to a peak flux density 
of around 70 mJy, which indeed corresponds 
to slightly less than 7o of our reported noise. 
The main result from these source counts is 
that we find evidence for flattening of the aver- 
age spectral index of the source population at 
23 MHz. The spectral index required to match 
the 54 MHz source counts from LoLSS [26] (- 
0.35+0.1) is flatter than the spectral index re- 
quired to match the 1.4 GHz source counts from 
to LoLSS (-0.5 +0.1, ref. [26]). For refer- 
ence, also contains the source counts 
from de Gasperin et al. (2021) [25], but scaled 
to the canonical spectral index of -0.7 instead 
(e.g. BIJ, illustrating that the source popula- 
tion reflects a flatter spectral index. This is 
in agreement with the findings in van Weeren 
et al. (2014) [32], where similar flattening at 
low frequencies has been found. Furthermore, 
this is also consistent with our finding that the 
median spectral index of our > 200 mJy inte- 
grated flux sources is —0.57 + 0.12 measured 
between 144 MHz and 23 MHz, which is also 
flatter than the canonical value of -0.7 between 
144 MHz and 1.4 GHz (see Methods). This 
curvature is also visible in Supplementary Fig- 
ure 4, where the majority of the nine brightest 
sources are shown to host curved spectra be- 
tween 1.4 GHz and 23 MHz, and none of these 
nine sources show inverted spectra. The curva- 
ture can be caused via absorption effects, which 
flattens the radio spectra at lower frequencies 
BIB], or via spectral aging, which causes steep- 
ening at higher frequencies [33]. 

Given that 98% of our detected sources are 
too faint to be detected by other decametric 
surveys, we compare our results to the 38 MHz 
8C survey, see Our results yield a 
significant improvement in both resolution and 
sensitivity over this survey. About 76% of the 
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Figure 4: Euclidean normalized source counts in the surveyed area (n = 2816, error bars corre- 
spond to expected value + standard error). Included are the source counts as constructed in this 
work, from Van Weeren et al. (2014) [32], from Mandal et al. (2021) [62], and from De Gasperin 
et al. (2021) [25], scaled with both a a = —0.7 and a = —0.35 power law. The corrected source 


counts are reported in Supplementary Table 1. 


sources cataloged by our LOFAR decameter ob- 
servations were undetected by the 8C survey, 
but are revealed in the LOFAR decameter ob- 
servations due to the significantly higher sensi- 
tivity. In addition, several cataloged individual 
8C sources are shown to be multiple sources, 
which are blended in 8C due to its worse reso- 
lution (4 arcmin). Of the 334 sources from 8C 
that are within the sky coverage of our obser- 
vations, 323 are detected by LOFAR and 11 are 
absent. Eight of the 11 undetected 8C sources 
are probably spuriously detected by 8C, given 
their low signal-to-noise ratio in the 8C cata- 
log. The other three sources are visible in the fi- 
nal image but were rejected in the catalog since 
these sources have a signal-to-noise ratio below 
seven. 


To exploit the first sensitive decameter 
23 MHz source catalog, we searched for sources 
having extremely steep spectra. We did this 
by cross-matching our data with the 151 MHz 
TGSS catalog and LoTSS. The TGSS survey 
covers the entire 305 square degrees area and 
achieves a similar depth to our 23 MHz catalog 
for compact radio sources with a typical AGN 
spectral index of —0.7. The search reveals two 


extended sources, not detected in TGSS, but 
visible in LoTSS, see [Figure 5| These sources 
also have the steepest spectral index among all 
the decameter sources detected in this work. 


The first source is located near the low-mass 
(Mso = 2.4 x 10!4Mo, ref. [34]) galaxy clus- 
ter Abell 553 (z = 0.067). Taking the 23 MHz 
flux density and the 144 MHz flux density from 
LoTSS, we compute an integrated spectral in- 
dex of Q23-144MHz = —1.82 + 0.13. For this 
calculation (and the others ones reported be- 
low), the 144 MHz flux density was extracted 
from a map tapered to the same resolution as 
the 23 MHz map. The steep-spectrum emission 
is located to the south and east of the brightest 
cluster galaxy (BCG). The source has a largest 
physical size of about 250 kpc. The LoTSS im- 
age shows an irregular morphology, with the 
brightest emission to the south of the BCG. 
A separate compact radio source is associated 
with the BCG itself. 


The second extended steep-spectrum source 
is located near the galaxy LEDA 2308099 (z = 
0.160, ref. [35]). It has a spectral index of 
Q23—-144 MHz = —1.43 + 0.13. The LoTSS im- 
age reveals a complex horseshoe-shaped source, 
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Figure 5: Images of four re-energised fossil plasma sources, detected at 23 MHz. These sources are 
associated with Abell 553 (top row), LEDA 2308099 (second row), Abell565 (third row) and Abell 566 
(bottom row). The left panel of each row contains a cutout made from LoTSS [24] data, with 23 MHz 
contours overlaid at levels of [8, 16,32] x o, with o = 14, 16, and 17 mJy beam’ for the top, second 
and third row respectively (corresponding with the local noise). The center panel contains a cutout 
of the 23 MHz map presented in this work, with TGSS [28] 151 MHz contours at [1,2,3] x o, with 
o =5mJy beam‘. The right panel contains a PanSTARRS gri band / Legacy image, overlaid with 
the same contours as the left panel. Due to the brightness of the source, the bottom left image has 23 
MHz contours overlaid at levels of [32,128,512] x c, with o = 18mJy beam™', and the bottom center 
image has TGSS 151 MHz contours at [4, 16, 64] x g- 


with the emission peaking to the west of 
LEDA 2308099. The radio emission has an 
extent of about 430kpc. A faint compact 
radio source is associated with the galaxy. 
LEDA 2308099 has no association with a known 
galaxy cluster or group. 

We also inspected resolved decameter sources 
with slightly less steep spectra. Most of these 
are classical radio galaxies. However, two of 
them have more peculiar morphologies. The 
first is located near the center of Abell 565, a 
low-mass cluster (M509 = 1.6 x 10!4Mo) at z = 
0.105 [86]. For this 320 kpc source, we measure 
a spectral index of @23—144 MHz = — 1.28 + 0.13. 
The source was previously studied by Mandal 
et al (2020) [87]. Higher frequency flux density 
measurements show that the spectrum steepens 
significantly to 325-1400 MHz ~ —2.05 + 0.08. 
The second source is associated with the low- 
mass [M500 = 3.0x10!4Mo, z = 0.098;[34] clus- 
ter Abell566. This complex 330 kpc extended 
source has a steep spectrum with a23_144 MHz = 
—1.26+ 0.12. It was also reported by Mandal 
et al. (2020) [37] and seems to be associated 
with the BCG. 

The irregular morphology, steep (curved) 
spectra, and presence of nearby radio AGN, but 
lack of clear direct association in the majority of 
cases, strongly suggest that the four sources dis- 
cussed above trace revived AGN fossil plasma 
38]. Revived fossil plasma (also known in the 
literature as re-energised fossil plasma or radio 
phoenices) is predicted to be produced by the 
re-acceleration of plasma from older AGN out- 
bursts. Simulations show that shock waves, for 
example from galaxy cluster mergers, can cause 
adiabatic compression of the old radio plasma 
[9] [12]. This produces radio sources with com- 
plex morphologies and steep, curved spectra. 
Very few such sources have so far been iden- 
tified [e.g., B9, [40], as their steep spectra make 
them faint at higher radio frequencies. Another 
expectation is that these sources are typically 
observed near cluster centers since the old radio 
plasma can be confined longer by the high den- 
sity of the intracluster medium. Furthermore, 
fossil sources should also be common in low- 
mass clusters, in contrast to larger megaparsec 
scale diffuse cluster radio sources [40]. Both of 
these expectations agree with our finding that 
three of our four detected sources are located 
near the centers of low-mass clusters. The dis- 
covery of these sources, with two of them being 


the steepest 23 MHz-detected source in the en- 
tire surveyed area, implies that fossil sources 
could form the dominant a < —1.5 popula- 
tion in the decameter band at our sensitivity 
limit. It is thus expected that many more can 
be found when more area is surveyed. 

From the results discussed in this work we 
conclude that sensitive sub-arcminute ground- 
based surveys of the decameter sky at fre- 
quencies above ~ 15 MHz are possible. Given 
the large field of view at these frequencies, 
such decametric surveys can be performed in 
a significantly shorter time than required for 
surveys at higher radio frequencies. More- 
over, ground-based observations are an attrac- 
tive option given that interferometers in space 
are costly. Finally, by optimizing a telescope 
for direction-dependent ionospheric calibration, 
it might be possible to push the limits to 
even lower frequencies than what is currently 
achievable with LOFAR. Previous work [41] has 
shown that in favorable conditions, the plasma 
cutoff frequency of the ionosphere can drop be- 
low 5 MHz. Given that the observations dis- 
cussed in this work were obtained under aver- 
age ionospheric activity, such an optimized tele- 
scope could perhaps operate at frequencies as 
low as 5-10 MHz during years around the solar 
minimum, when the ionospheric perturbations 
are the smallest. 


Methods 


Observation and initial data reduction 


For this work, an observation with LOFAR’s 
Low Band Antenna (LBA) in the 16-30 MHz 
frequency range was used. The observation 
was performed after midnight, as the impact 
of Radio Frequency Interference (RFI) at de- 
cameter wavelengths is significantly more se- 
vere during daytime and early evening than af- 
ter midnight (see [Figure 3). We used a 5 hour 
dual-beam observation, with beams centered 
on RA:06h20m15s, Dec.: +66°23m20s and on 
RA: 06h17m19s, Dec.: +54°20m49s; starting at 
00:00 UTC, 11 November 2021. Only the 
Dutch stations of LOFAR were used. With 
a third simultaneous beam, we observed the 
primary calibrator, 3C 196, which is located 
23°and 19°respectively from the center of the 
target beams. 


The data were reduced in a multi-tiered ap- 
proach. Firstly, we excised RFI from both our 
primary calibrator and target observation us- 
ing AOFlagger [42]. Next, the signals from 
the bright off-axis sources Cygnus A and Cas- 
siopeia A were removed from both the primary 
calibrator and the target observations [43]. We 
then applied a pre-determined bandpass and 
polarization alignment solution to the primary 
calibrator, after which we performed a phase 
and amplitude calibration to remove the ef- 
fects of the ionosphere and temporal gains. 
High-quality models of the primary calibrators 
sources at low frequencies are available for this 
purpose [44]. From the amplitude corrections, 
we determined the “S4-index” [45], which is 
a representation of the severity of the iono- 
spheric scintillations. The $4-index is defined 
as the standard deviation of the normalized sig- 
nal intensity in the ionosphere I(t)/ (I), and 
higher values correspond with a more active 
ionosphere. In order to verify that our obser- 
vation represents a normal ionosphere, we ana- 
lyzed observations of our primary calibrator on 
41 other nights with different ionospheric con- 
ditions between October and December 2021. 
The $4-index of these observations was between 
0.09 and 0.5, with a median of 0.20 and a stan- 
dard deviation of 0.10. The $4-index of the ob- 
servation presented in this work was 0.18, which 
is comparable to the median value of the set of 
calibrator observations. 

Subsequently, we applied the same bandpass 
and polarization alignment solutions to the tar- 
get fields, as well as the phase and amplitude 
solutions obtained from the primary calibrator 
3C 196. This step effectively takes out the in- 
strumental effects, leaving only ionospheric per- 
turbations. These perturbations have a clear 
frequency-dependent functional form. 

Next, we self-calibrate the full field of view 
of each target pointing separately, using the 
74MHz VLSSr [46] as a starting model, cor- 
recting for the average TEC effect across the 
full image and differential Faraday Rotation in 
the ionosphere. Differential Faraday rotation is 
taken out by first transforming the visibilities 
to a circular basis, where this effect becomes a 
phase difference between the RR and LL cor- 
relations. We then solve for this phase dif- 
ference using a point source model, employing 
DP3 [47]. This is followed by an antenna-based 
phase-only self-calibration cycle starting from 
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the VLSSr model, constraining the phase cor- 
rections to be smooth along the frequency axis. 
This allowed us to make a preliminary image of 
the target field with WSClean [48], facilitating 
the identification of the brightest objects in the 
field of view. 


Direction-Dependent calibration 


As the ionosphere differs significantly across 
the field, most sources in the image with the 
direction-independent calibration are still sig- 
nificantly distorted. For this reason, the field 
is broken up in facets, which are subsequently 
self-calibrated, with a technique similar to de 
Gasperin et al. (2021) [25] and Van Weeren 
et al. (2016) [27]. For the self-calibration we 
make use of the facetselfcal script described 
in Van Weeren et al. (2021) [49], solving for dif- 
ferential Faraday rotation and antenna-based 
phases with DP3 that are constrained to be 
smooth along the frequency axis. As a start- 
ing model, we use TGSS-ADR1 [28]. This could 
also be LoLSS or LoTSS-based model, however, 
these were not yet available for the entire region 
of sky covered by our observation. Initially, this 
calibration step is only possible for facets con- 
taining the brightest few (> 10 Jy) sources, as 
fainter sources are severely perturbed by un- 
corrected ionospheric effects from the brightest 
nearby sources. Using DDFACET [50], we are 
able to create an image and model, correct- 
ing each direction during imaging. This im- 
proved model allows us to break up the field in 
more facets (about a factor 2 times more), as 
the ionospheric effects from bright sources are 
now better corrected. The increased number of 
facets that can be calibrated results in an im- 
age that is sharp across a larger fraction of the 
field of view, allowing us to iteratively increase 
the number of facets that can be calibrated. 
This cycle is repeated for in total 4 direction- 
dependent calibration cycles, giving us a total 
of 25 facets across the field of view of a single 
target pointing. 

For each target pointing, the image was re- 
stored with a circular 45 arcsecond beam, which 
is close to the native resolution. From this im- 
age, we extracted sources with PyBDSF [29] into 
a preliminary catalog. Due to residual iono- 
spheric effects (particularly caused by residual 
scintillation and direction-dependent errors), 
unresolved sources are not well approximated 


by a Gaussian profile, which means that the fit- 
ted fluxes could be inaccurate. For this reason, 
circular apertures centered on the positions ob- 
tained with PyBDSF, with a radius of 2 arcmin- 
utes are used to measure the flux for each source 
instead. If sources are located within 4 ar- 
cminutes of each other, the corresponding aper- 
tures are truncated by the bisector of the line 
segment connecting the sources, ensuring that 
each source contains the flux corresponding to 
the sky area closest to the source. Within each 
aperture, the flux density is computed by sum- 
ming the intensity inside each aperture and di- 
viding it by the area of the restoring beam. 
Only sources with a ratio of peak flux den- 
sity of 5 times the local RMS, an integrated 
flux density noise ratio above 7 and a maxi- 
mum distance to the pointing center of 7 de- 
grees were considered. As this approach dif- 
fers from the more conventional approach for 
creating radio catalogs [e.g. [24] [26], we com- 
pare the flux densities obtained in this work 
with the raw fluxes created with PyBDSF. In 
addition, spurious sources located on calibra- 
tion artifacts near bright sources were manu- 
ally excised. We see in Supplementary Figure 
1 that PyBDSF somewhat underestimates the 
flux densities of sources compared to aperture 
flux densities, however this underestimation is 
less pronounced for brighter sources, as the cal- 
ibration solutions of the facets are weighted to- 
wards these brighter sources. This is expected, 
as residual artifacts are not effectively fitted in 
PyBDSF but are included in apertures. 

The aforementioned residual ionospheric ef- 
fects may cause blurring of sources, which 
causes a reduction in effective resolution. This 
results in a reduction of the peak flux den- 
sity compared to the integrated flux density, as 
the blurring causes the full width at half max- 
imum (FWHM) to be increased. Supplemen- 
tary Figure 2 shows the ratio of integrated flux 
density over peak flux density, where the ratio 
of integrated flux density over peak flux den- 
sity corresponds to the broadening caused by 
residual ionospheric effects. For faint sources 
with a peak over local RMS ratio of around 
~ 5, peak flux densities are reduced by a factor 
~ 2.45 compared to the integrated flux densi- 
ties, which means that the ionospheric blurring 
is around ,/(2.45- 45)? — 452 = 100”, while 
for bright sources, the ratio of peak flux den- 
sity over local RMS is around ~ 1.67, which 
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means that the ionospheric blurring is around 
y (1.67 - 45)? — 45? = 60”. 


In order to verify the astrometric accuracy 
of our observation, we compared our catalog 
with the TGSS catalog at 151 MHz, which has 
a comparable resolution to the data presented 
in this paper. In Supplementary Figure 3, it 
can be seen that the astrometric offsets are sig- 
nificantly smaller than the resolution of the LO- 
FAR data, which indicates that the astrometry 
is reliable. We also compare the flux densities 
measured by TGSS with the flux densities mea- 
sured at 23 MHz in Supplementary Figure 4. 
In Supplementary Figure 5, we show the flux 
density of 9 sources with the highest flux den- 
sity detected with LOFAR below 30 MHz, com- 
pared to previous source catalogs, to show that 
the flux density is in line with expected val- 
ues given higher frequency data. This agrees 
with de Gasperin et al. (2021) [25], who also 
find that the LOFAR LBA flux-scale can be re- 
liably transferred from the calibrator beam to 
the target beam. 


Data Availability 


The catalog produced in this work is available in the Supplementary Data. The calibrated mosaic 
image is available upon reasonable request to the authors. 


Code Availability 


The following codes are used in this work, and readily available: 
e https://github.com/saopicc/DDFacet (DDFacet) 
e https://git.astron.nl/RD/DP3 (DP3) 
e https://gitlab.com/aroffringa/wsclean (WSClean) 
e https://github.com/rvweeren/lofar_facet_selfcal (facetselfcal) 
e https://github.com/lofar-astron/PyBDSF (PyBDSF) 


Other code used in this work is available upon request to the authors. 
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